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Abstract 

Fires and rainfall are major mechanisms that regulate woody and grassy biomasses 
in savanna ecosystems. Conditions of long-lasting coexistence of trees and grasses 
have been mainly studied using continuous-time modelling of tree-grass competition. 

In these frameworks, fire is a time-continuous forcing while the relationship between 
woody plant size and fire-sensitivity is not systematically considered. In this paper, 
we propose a new mathematical framework to model tree-grass interaction that takes 
into account both the discrete nature of fire occurrence and size-dependent fire sen¬ 
sitivity (via two classes of woody plants). We carry out a qualitative analysis that 
highlights ecological thresholds and bifurcations parameters that shape the dynamics 
of the savanna-like systems within the main ecological zones. Moreover, through a 
qualitative analysis, we show that the impulsive modelling of fire occurrences leads 
to more diverse behaviors and a more realistic array of solutions than the analogous 
time-continuous fire models. Numerical simulations are provided to illustrate the theo¬ 
retical results and to support a discussion about the bifurcation parameters and future 
developments. 

key words: Asymmetric competition - Savanna - Fire - Impulsive differential equation - 
Qualitative analysis - Nonstandard finite difference scheme 

1 Introduction 

Savannas are ecosystems with fairly continuous grass cover and variable woody cover (Mau- 
rin et ah (2014) |31])- However, savanna-like ecosystems are diverse and cover extensive 
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areas throughout the tropics. Explanations found in the literature about the possible long- 
lasting coexistence of woody and grassy vegetation components therefore relate to diverse 
factors and processes depending on the location and the ecological context (Baudena et ah 
(2014) [n]). Several studies have pointed towards the role of stable ecological factors e.g. 
climate, in shaping the tree to grass ratio along large-scale gradients of rainfall or soil fer¬ 
tility (Sankaran et ah (2005) |10], (2008) [39]). Other studies have rather emphasized the 
reaction of vegetation to recurrent disturbances such as herbivory or Ere (Langevelde et ah 
(2003) [21], D’Odorico et ah (2006) [20], Sankaran et ah (2008) [20], Smit et ah (2010) |4T] . 
Favier et al. (2012) [21] and references therein). Those two points of view are not mutually- 
exclusive since both environmental control and disturbances may co-occur in a given area 
and along ecological gradients, although their relative importance generally varies among 
ecosystems. Bond et al. (2003) [H] proposed the name of climate-dependent for ecosystems 
which physiognomies are highly dependent on climatic conditions (rainfall, soil moisture) 
versus disturbance-dependent for ecosystems which dynamics is strongly dependent on hres 
or herbivores. 

Several models using a system of ordinary differential equations (ODES) have been pro¬ 
posed to depict and understand the dynamics of woody and herbaceous components in 
savanna-like vegetation. A Erst attempt (Walker et al. (1981) ^U]) was orientated towards 
semiarid Ereless savannas and analyzed the efiect of herbivory and drought on the balance 
between woody and herbaceous biomass. This model refers to ecosystems immune to Ere due 
to insufEcient annual rainfall and grass production. Indeed, Eres in savanna-like ecosystems 
mostly rely on herbaceous biomass that has dried up during the dry season. As long as 
rainfall is sufEcient, Eres impact seedlings and saplings within the Eame zone and thus let 
grasses indirectly inhibit tree establishment. 

More recently, several attempts have been made (see Langevelde et al. (2003) [3T] . 
Accatino et al. (2010) [2], De Michele et al. (2011) [19], Tchuinte et al. (2014) [4^, Yatat et 
al. (2014) [22]) to model the dynamics of savannas, taking into account Eres as continuous 
events, on the basis of the initial framework of Tilman (1994) [46] that used coupled ODES 
to model the competitive interactions between two kinds of plants. 

However, it is questionable to model Ere as a continuous forcing that continuously removes 
fractions of Ere sensitive biomass. Indeed, several months and even years can past between 
two successive Eres, such that Ere may be considered as an instantaneous perturbation of 
the savanna ecosystem. Several recent papers have proposed to model Eres as stochastic 
events while keeping the continuous-time difierential equation framework (Baudena et al. 
(2010) [TU], Beckage et al. (2011) [12]) or using time discrete matrix models (Accatino & 
De Michele (2013) [2]). But in all those examples. Ere characteristics remain mainly a linear 
function of grass biomass which is not satisfactory. Indeed, it is well known that at low 
grass biomass there is no Eres while above a sufEcient grass biomass, Eres intensity increases 
rapidly before reaching a saturation. This particular feature of Eres in savanna vegetation 
cannot be modeled by a linear function. Another drawback of the aforementioned recent 
stochastic models (Baudena et al. (2010) [10], Beckage et al. (2011) [12]) is that they barely 
lend themselves to analytical (qualitative) approaches. 

In this paper, we therefore present a model that difiers from most published models and 
extends the work of Yatat et al. (2014) [22] by modelling discrete Ere occurrences. We 
consider a tree-grass compartmental model. We set one compartment for grass biomass and 
two for trees, namely Ere-sensitive individuals having most of their buds within the Eame 
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zone (like seedlings, saplings, shrubs) and non-sensitive mature trees having at least their 
upperparts above the flame zone. We therefore develop a system of three coupled non-linear 
impulsive differential equations (IDES), one equation per vegetation compartment, that 
describes savanna dynamics. In addition, we model Are intensity (and the corresponding 
impact on sensitive woody plants) as an increasing nonlinear and bounded function of grass 
biomass. Finally, Are occurrences are modeled as pulse-like perturbations. 

In order to assess and illustrate the dynamics of some ecological formations through 
our mathematical model, we distinguish in our numerical computations three climatic zones 
having distinct characteristics in terms of biomass production. These biomass production 
zones, loosely relate to Africa and are indexed by their carrying capacity for grass and woody 
biomasses and by fires return times. Semi-arid areas have a mean annual rainfall that varies 
between 300 mm.yr~^ and 650 mm.yr~^, and fires, if any have return intervals relatively 
long, says less than one Are event every ten years. Mesic savannas have a mean annual 
rainfall comprised between 650 mm.yr~^ and 1100 mm.yr~^ and Are return time is in order 
of four or flve years, sometimes less. Finally, we consider a humid tropical area with a mean 
annual rainfall between 1100 mm.j/r“^-1800 mm.yr~^ and where one can have a Are return 
time from less than one year to two-three years. 

Although impulsive diAerential equations appear highly relevant to model vegetation dy¬ 
namics in Are prone savannas, they are also difhcult in terms of analytical treatments. This 
may explain why this framework has still remained scarcely used for modelling the dynam¬ 
ics of Are prone savannas. Our model aims to acknowledge three major phenomena: the 
periodic occurrence of Are events, the Are-mediated, non-linear negative feedback of grasses 
onto sensitive trees and the negative and/or positive feed-back, depending on location, of 
insensitive trees on grasses. We therefore explicitly model the occurrence of Ares in savanna 
ecosystems and the asymmetric nature of tree-grass competitive interactions in savannas. 

The full impulse Are model of asymmetric tree-grass competition (IFAC) is formulated 
in Section [2j In Section |3] we reach qualitative analytical results for IFAC through which 
we highlight some meaningful ecological thresholds that summarize savanna dynamics under 
impulsive Ares. We present a nonstandard numerical scheme for the IFAC model in Section 
m together with the IFAC parameters ranges. In section [5] we present numerical simulations 
done in the three ecological biomass production areas of the African continent. Moreover, 
Section [5] also deals with the discussion of our results. 

2 The impulse fire model of asymmetric tree-grass com¬ 
petition (IFAC) formulation 

As we have mentioned before, we consider vegetation as composed of three classes, i.e. the 
class of sensitive tree biomass (T^), the class of non-sensitive tree biomass (T/vs) and the 
class of grass biomass (G). We model the Are intensity by a non-linear increasing function 
of grass biomass w{G). To built up our model, we consider the following assumptions where 
(A4), (A5), (A6) and (A7) are already described in Yatat et ah (2014) [52] . 

(Al) A carrying capacity Kt for tree biomass (in tons per hectare, t.ha ^). 

(A2) A carrying capacity Kq for grass biomass (in tons per hectare, t.ha ^). 
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(A3) Fire events occur period.cally, ..e. every r-tune, where r = - and / denotes the fire 
frequency. 

(A4) Fire only impacts grass and sensitive Tree and, fire intensity is an increasing function 
of the grass biomass 


(A5) The Grass biomass has a direct, depressing effect on the Sensitive Trees. 

(A6) Non Sensitive Trees have a depressive or facilitation effect on grass biomass by shading. 

(A7) Sensitive tree biomass moves to non-sensitive tree biomass after an averaged time — 

Us 

(in years). 

The following diagram summarize the relationship between the three compartments 
75 



Figure 1: Compartmental diagram of the size structured tree-grass interaction model in 
impulse hres-prone savanna. 


The following parameters are considered throughout the paper: 

• Sensitive tree biomass is made up from non sensitive tree biomass (i.e. seed production 

and germination) with the rate 'Jns and from existing sensitive tree biomass 

(i.e. intrinsic growth) with the rate 75 {in yr~^). 

• Grass biomass is made up from existing grass biomass with the intrinsic growth rate 
7g {in yr-^). 

• /i 5 {in yr~^) is an additional death rate of sensitive tree biomass due to external 
disturbances such as human activities and herbivory. 

• yNS {'^^ y'^~^) is fhe natural death rate of non sensitive tree biomass. 
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• / is the fire frequency (in yr ^). 

• hg {in yr~^) is an additional death rate of grass biomass due to factors including 
human activities and herbivory. 

• — (in yr) is the average time that a sensitive tree takes to become non sensitive to 
Us 

hre. 


• ac expresses the asymmetric competition exerted by grasses on sensitive trees (shading 
and competition for nutrients, in ha.t~^.yr~^). 

• o'Ns expresses the asymmetric interaction (competition and/or facilitation) of non sen¬ 
sitive trees on grasses (shading and competition for nutrients, in ha.t~^.yr~^). 

• 7^5 is the proportion of sensitive tree biomass that is consumed by hre. 

• rjG is the proportion of grass biomass that is consumed by hre. 


Based on these ecological premises, and taking into account the ehect of hre as pulse 
phenomena, we propose a model for the savanna vegetation dynamics which is and extension 
of the model formulated in Yatat et al. (2014) [52] . 

The IFAC is given by 


dt 

dT^s 

dt 

dt 


with 


{isTs + InsTns) ( 1-1 — TsifJ'S + + (^gG), 


K. 


^^sTs — f'NsTNS, 
G 


t ^ tk (1) 


7g ^ ~ {^nsTns + t‘'G)G, 

Tsitt) = {l-ysw{G{t,)))Ts{h), 

TNs{ik) = T^sitk), t = tk tk+i = tk + T 

G{tt) = {l-yG)G{tk), 

Ts{0) = Tso > 0, Tns{0) = TnSo ^ 0 ) ^(0) = GO > 0. 


( 2 ) 

(3) 


For this impulsive hre model, the hre intensity function tc is a continuous and positive 
function of grass biomass which is bounded above by unity. As in Tchuinte et al. (2014) [45] 
and Yatat et al. (2014) [52], one can choose a generic sigmoidal function (see also Staver et 
al. (2011) [13]). A typical choice could be 


m;(G) 




(4) 


where Gq = is the value of grass biomass at which hre intensity reaches its half saturation 
{qq in tons per hectare, t.ha~^) and a G N*. 

The feasible region for system ([I]) — ([2]) is the set kl dehned by 


kl = {{Ts-,T^s-,G)&^1\Q<Ts + T^s<Kt,Q<G<Kg}. (5) 
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3 Mathematical Analysis and Ecological Interpretation 
of Thresholds 

3.1 Existence of solution 

The right-hand side of system ([ID-([2]) is locally lipschitz continuous on f]. Thus, using a 
classic existence Theorem (Theorem 1.1 P. 3 in Bainov and Simeonov (1995) [8]), system 
©- 0-0 has a unique solution on 

3.2 Trivial and semitrivial solutions 

It is obvious that system ([I]) — ([2]) has always a desert equilibrium 

Eo = (0,0,0). 

3.2.1 The positive grassland periodic solution: existence and local stability 


Let us consider the following thresholds: 

K = —, f^G> 0, (6) 

f^G 

and 

nO = / - l) ) /^G > 0, , . 

1 ( 1 - Vg) exp (jgt) , /iG = 0. ^ ^ 

Assume that TZq > 1, we have the following result (see also Dai et ah (2012) [18]) 

Lemma 1 When Pq > 1, System 0-0 has a positive grassland periodic solution Eg = 
(0; 0; G*(t)), where 


G*(t) 


Kg {pI - 1) 

(p^ - 1) + ’ 


Pg > 0 , 


nr <t < (n + 1)t, 


Pg — 0 . 


( 8 ) 


Remark 1 (Thresholds interpretation) • TZq is the average amount of biomass pro¬ 
duced per unit of grass biomass during its whole lifespan in absence of fires and de¬ 
pressing effect from non sensitive trees but subject to additional mortality caused by 
human activities or by herbivory. 


• p'q embodies the residual amount of grass biomass at any time-period that fires occur, 
from the grass biomass produced per unit of grass biomass. 


Remark 2 An obvious computation leads to: 


< 1 
Pg > 1 
P°G>1 



Pg < 

Pg < 1 - exp(-pG(^G 

Pg < 1 - exp(-7Gr), 


l)r), Pg > 0, 

Pg = 0. 


(9) 


6 






Let us state the following Lemma which will be helpful for the sequel 


Lemma 2 When > 1, 


'•(n+l)T 


Gint — 


T 


G*{s)ds = hi(p^) = -— (ln(l - pc) + (7 g - > 0. 

r 7 g r 7 g 


( 10 ) 


Now we turn to look for local stability of that previous positive and periodic grassland 
solution. For that purpose, we will use the small perturbation technique and the Floquet 
theory, i.e. we will hud conditions under which all the Floquet multipliers of the positive 
and periodic grassland solution have their absolute value less than unity or equal to unity 
(D’Onofrio (2002) [2l], Chen et al. (2009) [I6]). 

Let 

,T 'ysh'NS + l^S'lNS 


I^Nsil^S + l^s) + I^NSO’cGint 


ISh^NS + ^^SlNS 


( 11 ) 


hvs(hs + (^s) + (ln(l - 7g) + (7g - hG)r) 

T IG 


^ = 75r ( 1 - - ) , 


where 


and 


7^ = 


_7s_ 

hs + + h'NS + CrcGint 


fS — (hs + l^S + O'GGint) (l — Tic) ■ 

Moreover, let Ai, A 2 be the roots of 

p(A) = \^ - AX+ B 


and 


Pt = max{(l - risw{G*{T)))e^\ e^"}. (12) 

The following result holds for System ([I]) — ([2]). 

Lemma 3 (Local stability of the grassland periodic solution Eg) 

Suppose that the grassland periodic solution (Eg) exists i.e. 

TZq > 1 and Pq > 1. 

Moreover, 

• ifR-G < 1 then, Eg is locally asymptotically stable, 

• if (TZg > 1 and pr < 1) then. Eg is locally asymptotically stable, 

• if (JIq > 1 and pt = 1) then. Eg is locally stable. 

• if (TZq > 1 and pt > 1) then. Eg is unstable. 
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Proof 1 (See Appendix A.) 


Remark 3 (Thresholds interpretation) In the sequel, we provide approximative thresh¬ 
olds interpretation in order to favor an intuitive ecological comprehension of our results with 
respect to these thresholds. 

• IZq is the sum of the average amount of biomass produced by a sensitive/young plant 
competing with grass between two successive fires, and the average amount of biomass 
produced by a mature plant multiplied by the proportion of young plants which reach 
the mature stage. 

• pt embodies both the residual of the reduction of trees biomass due to periodic fires 
events and the reduction of sensitive tree biomass due to competition with grass biomass. 
Moreover, since IZ < after a direct computation one has 

ac < —{'Is — ih'S + + h'Ns)) Pt > 1- (13) 

Therefore, it clearly appears following relation / f73j) that the grass vs. sensitive tree 
competition parameter ac is a bifurcation parameter for the IFAC model that embodies 
the stability/instability of the grassland periodic solution. 

3.2.2 The positive forest equilibrium: existence and local stability 


Let 


= 


Tspvs + iNS^^S 


Tns{ts + ^s) 

The following result follows from Proposition 1 in Yatat et al. (2014) 


(14) 


Lemma 4 IfIZf > 1 then system ([I]) — ([2]) has a positive forest equilibrium Et = {Tg^T^s]^), 
where 


Ts = 
Tns = 


KtTns 
Pns + 


1 - 


^^-fs = 


Ktojs 


n% 

1 - 


TTf 


(15) 


Pns Pns + ^s 

As previously we are checking for local stability of the positive forest equilibrium. Let 

^T = - ^ , 

Pg + o'nsTns, 

JG _ 


Pt = (1 - Pg) exp (^1 - ^ J r 

Using the same approach as in the proof of Lemma [31 we derive the following result. 

Lemma 5 (Local stability of the forest equilibrium Et) 

Suppose that the forest equilibrium (Et) exists, i.e. IZf > 1. 

Moreover, 


( 16 ) 


if^Z/f ^ 1 then, E/ is locally asymptotically stable. 








• if (JZ^ > 1 and < 1) then, Et is locally asymptotically stable, 

• if (JZ^ > 1 and = 1) then, is locally stable . 

• if (JZ^ > 1 and p^ > 1) then, Et is unstable. 


Remark 4 (Thresholds interpretation) As we mentioned before, we provide approxima¬ 
tive thresholds interpretation in order to favor an intuitive ecological comprehension of our 
results with respect to these thresholds. 

• TZf is the sum of the average amount of biomass produced by a sensitive/young plant, 
without fires and competition from grass, and the average amount of biomass produced 
by a mature plant multiplied by the proportion of young plants which reach the mature 
stage. 

We may note here that this threshold only depends on the parameters ruling the dy¬ 
namics of the woody biomass. 

• TZt is the average biomass produced by a unit of grass biomass during its whole lifespan 
free of fires while experiencing competition from non-sensitive trees. 


• Pt embodies both the residual grass biomass after a fire event and the residual biomass 
of the depression of grass biomass due to competition from non-sensitive trees. 

One should also note that, when TZ/ > I, one has 


Pr < 1 Pg > 1 


exp ( 7 g ( 1 - ^ ) r 


r < 


ln(l - pg ) 




(17) 


1 f , ln(l - Pg) , ^ 

=— 7g - Pg +- 1 < (^NS- 

NS \ ^ 


Therefore we can deduce three major observations 


(i) Firstly, if the fires period r is small (i.e. the fires freguency f is high) then 
1- 7-7 -;— r —^ is Small and one can have pS < 1 for small values of 

-P(-i') 

Vg- 

(a) Secondly, if the fires period r is large (i.e. the fires freguency f is small) then 


where 0 < e << 1. Thus, e can be sufficiently small (such that Pt <l), to have 
a large destruction of grass biomass (more than 99%). This is ecologically not 
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possible (according to the fact that a part of the grass biomass like roots and even 
the bottom of the tufts, cannot burn). Therefore, in case of large fires period, 
having < 1 may likely correspond to a decrease of (i.e. an increase of 
(^Ns)- 


(in) It is easily deduced from relation l[T7\ ) that the non sensitive tree vs. grass competi¬ 
tion/facilitation parameter a^s the fire return time r are bifurcation parame¬ 
ters for the IFAC model that embody stability/instability of the forest eguilibrium. 


Remark 5 A direct comparison leads to: 

1 . 

2. Tic <l k‘G> 0. 


3.2.3 Global stability of trivial equilibrium (desert) and semi-trivial solutions 
(grassland periodic solution and the forest equilibrium) 


Here, we state a result concerning the global stability of the desert equilibrium, the forest 

equilibrium and a result concerning the global stability of the grassland periodic solution. 

Using the thresholds dehned in ([S]), ([7]) and ([HD, we have the following 

Theorem 1 • Case 1: nc > 0. 

1. IfHq^ < 1 and TZq < 1 then the desert eguilibrium Eq is globally asymptotically stable 
(GAS). 

2. IfHq- > 1 and TZq < 1 then the forest eguilibrium Et = (Ts,T/v5,0), where {Ts,T]ys) 
are given in m. is globally asymptotically stable. 

3. ifn% ^ 1, Hq ^ 1 and Pc ^ 1 then the desert egu%l%br%um Eq %s globally asymptotically 
stable. 

I ifK ^ 1, He ^ 1 and Pc ^ 1 then the grassland periodic solution Ec (0, 0, G (t)), 
where G*(t) is given by is globally asymptotically stable. 

• Case 2: pc = 

(i) IflZf < 1 and < 1 then the desert eguilibrium Eq is globally asymptotically stable. 

(a) If Hff < 1 and p^ > 1 then the grassland periodic solution Ec = (0,0,G*(t)), where 
G*(t) is given by is globally asymptotically stable. 

(Hi) IfHj- > 1 and p^ < 1 then the forest eguilibrium Et = (T 5 , Tats, 0), where {Ts,Tits) 
are given in m, is globally asymptotically stable. 

Proof 2 (See Appendix B.) 
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3.3 Existence of a positive and periodic tree-grass solution 

Now, we reach the position to hnd at least one non-trivial positive and periodic solution of 
system (HD-©. We will use the approach developed by Gaines and Mahwin 1977 [27] . 
Before we give the main result of this section, we recall useful inequalities according to the 
thresholds dehned in ([ 6 ]) and ([7]) 

1. If TZ^q < 1 then < 1. 

2. > 1 ^ (7g - Pg) + ~ > 0. 

The following result holds for system © — ©. 

Theorem 2 (Existence of a positive and periodic savanna solution) 

• Case 1: pc > 0. 

IfE-Q > 1 and p^ > 1 then system © — © has at least one positive t- periodic solution. 

• Case 2: pc = 0. 

If Pq > 1 then system © — © has at least one positive r-periodic solution. 

Proof 3 (See Appendix C.) 

Remark 6 One should note that Theorem provides only sufficient conditions to ensure 
existence of at least one positive r-periodic solution of system In other words, 

existence of savanna solution relies on sufficient grass biomass production. Thus, from an 
ecological point of view assumptions of Theorem \E can also be view as necessary conditions 
to have savanna solution. 


Remark 7 Note also, that the unigueness of the positive savanna solution is an open prob¬ 
lem. In the rest of the paper, we will assume that we only have one positive solution. 


As previously we can check the local asymptotic stability of the positive and periodic 
savanna solution Etg = iTs{t),Tj^s{t),G{t)). Dehning 


Ts{l) = x{t)+Tj{t), 
Tns{ 1) = y{l) + TNs{l)i 
Git) = zit)+Git), 

where a;(f),p(f) and zf) are small perturbations and satisfy 


xf) 


a;(0) 


y(t) = ®(«) !/( 0 ) 


z(t) 


2 ( 0 ) 


(18) 


(19) 


where $ is a fundamental matrix and satishes 


d^jt) 

dt 


,(i) 




a( 2 ) 

-PNS 


0 —o'NsGit) 



m ( 20 ) 
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with 


= 75 1^1 

= Ins f 1 - 

a(3) = I 1 


_ ‘2Ts{t) + TNs(t) \ _ Ins 
Kt j Kt 


Kt 

Tsit) + 2T7V5(^) 


TNs{t) — O'GGit) — fis — <^s, 


Ki 


- 


( 21 ) 


20 ( 1 ) ^ 
Ka / 


O'NsTNsi't) — he 


and *h( 0 ) = JdRS. Furthermore, the resetting impulsive condition of system (II])-(I 2 ]) becomes, 


/ x{nT+) \ / 1 - r]sw{G{T)) 0 -w'(G(r))Ts(r) \ / x{nT) \ 

y{nT+) = 0 1 0 yinr) . (22) 

\ z{nT+) / \ 0 0 1-7G / \ z{nT) / 

A monodromy matrix M of system ffTl) — l|2|), is: 

/ 1 - ?7sw(G(r)) 0 -w'(G(r))Ts(r) \ 

M = 0 1 0 <h(r), (23) 

V 0 0 l-yc J 

with ^ 

d)(f) = exp(^j^ DF{fs{s);fNs{s);G{s))dsy (24) 

Let dehned ptg such as 

Ptg = max(|A| : A G sp{M)), (25) 

where the matrix M is dehned in fl25]l . Therefore, following the Floquet theorem (D’Onofrio 
(2002) [2T], Chen et ah (2009) [I6j) we deduce the following results 


Lemma 6 • If Ptg < 1 then the positive r-periodic solution of system ([I]) — ([2]) is locally 

asymptotically stable. 


• If Ptg = 1 then the positive t- periodic solution of system ([I]) — ([2]) is locally stable. 

• If Ptg > 1 then the positive t- periodic solution of system ([I]) — ([2]) is unstable. 

Unfortunately, expressions fl23l) and fl2T)l don’t allow an explicit computation of the eigenval¬ 
ues of the monodromy matrix M and of the real pxG- Therefore, the stability of the positive 
r-periodic solution of system ([I]) — ([2]) will be conjectured through numerical computation 
of the threshold ptg- 


3.4 Summary Table of the IFAC model qualitative analysis 

Based on the previous studies, we deduce the summary Table of the qualitative analysis of 
the IFAC model. In Table 13.41 L stands for Locally Asymptotically Stable, NU stands for 
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Numerical Asymptotical Stability and the empty cell denotes either the instability/non exis¬ 
tence of the corresponding solution or that the result does not depend on the corresponding 
threshold. For reader’s convenience, we recall all thresholds dehned previously: 


11 % 


IS^NS + iNS^S 




—, when /ic > 0 


Pg 


= (1 - r]G) exp (pg - 1 ) r), 


^ , 

Pg + ctnsTns 


where Tns = W 

Pns + V 


Pt 


(1 - //g) exp ( 7 g ( 1 - ^ I r 




IsPnS + 


PNsiPs + (^s) + — (ln(l - Pg) + {ig - Pg)t) 

T 7 g 


and Pt is dehned in ([ 12 ]). 

In Table 1231 we implicitly assume that Rq > 1 and R% > 1 . 


(26) 


Thresholds 

Solutions 


n% 

Pg 

Pt 

nl 

Pt 

Et 

Eg 

Etg 

Case 

> 1 

> 1 

> 1 

> 1 

> 1 



NU 

1 





< 1 


L 

NU 

2 




< 1 



L 

NU 

3 



< 1 

> 1 

> 1 

L 


NU 

4 





< 1 

L 

L 

NU 

5 




< 1 


L 

L 

NU 

6 


< 1 




L 


NU 

7 

< 1 

> 1 


> 1 

> 1 

L 


NU 

8 





< 1 

L 

L 

NU 

9 




< 1 


L 

L 

NU 

10 


< 1 




L 


NU 

11 


Table 1: Summary table of the qualitative analysis of system ([I]) — (|2|). NU: LAS for Etc 
needs to be estimated numerically, using 


Remark 8 In Table we only summarizes the cases where TZq > 1 and 1Z% > 1, which 
are the interesting cases from the ecological point of view. Moreover, a direct computation 
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leads 


Pt = (1 - Vg) exp ^ j rj 

= {l-r]G)exp{{'yG-PG-(7NsfNs)r) 

= (1 - Vg) exp (pg('^g “ 1 )^) exp {-ctnsTnst) 


(27) 


Pt = Pc^wi-'^CTNsTNs)- 

Therefore, aeeording to ( [F7| ), since TZf >1 i.e. 1^5 > 0 then < p^. 

In the sequel, we provide some numerical simulations in order to illustrate our theoretical 
results. To achieve that goal, hrst we will provide a suitable nonstandard numerical scheme 
which will be helpful for the numerical approximation of the IFAC model solutions. It is 
well know that standard methods (such as Runge-Kutta or Euler methods) can sometimes 
present spurious behaviors which are not in adequacy with the system properties that they 
aim to approximate i.e., lead to negative solutions, exhibit numerical instabilities, or even 
converge to the wrong equilibrium for certain values of the time discretization or the model 
parameters (interested readers can also see Yatat et al. (2014) [52], Anguelov et ah (2012) 
[I], (2013) 0, (2014) [6| and Dumont et al. (2010) [22], (2012) [23] for motivations, details 
and explanations about nonstandard schemes). Secondly, we will focus on three ecological 
regions of the African continent that contrast in terms of biomass production conditions, 
namely a semiarid, a mesic and a humid tropical region, to discuss the IFAC outcomes with 
respect to published modelling results on savanna ecosystems (Baudena et al. (2014) [TT] . 
February et al. (2013) [25], Accatino et al. (2010) [3], Mordelet et al. (1995) [36], Moustakas 
et al. (2013) [37]). 


4 Nonstandard scheme, parameters ranges and ecolog¬ 
ical zones of biomass productions 

4.1 A nonstandard scheme for the IFAC model 


The nonstandard numerical scheme proposed in this section is adapted from the nonstandard 
scheme proposed for the COFAC model in Yatat et al. (2014) [52] . 

System ([ID is discretized as follows: 


^fc+l _ Qk 

f^G^h) 


= 7g 1 - 


Tlfc + l 

-‘-NS 


rpk 

-‘ns 


0 (/l) 


rj-^k + l 


rpk 

-‘s 


(j){h) 


Gk+1 


K, 


G 


- aNs^sG^^^ - PgG\ 


= ^sTs - pnsT^\^, 


(75 - {ps + 0Js))Tl + iNsT^Nf - ir'^siTs 


'yNS rjik rTik-\-l ( 'yNS rjik I /^k \ rrifc + l 

-i ATQ-I- I -^NS “I" I -^,9 


'yS rjik f rjik-\-l I rjik-)-!^^ 

“T -L i 


NS 




'S ’ 


( 28 ) 
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and the impulsive event ([ 2 ]) is discretized as follows: 


, rpk+l _ rpk+1 

^ -‘-NS+ — -‘-NS ^ 

. Tit' = (1 - >)s«.(G'=+'))rp‘, 

where the denominator functions 0 and 0 i read as 

Qh _ 

(p{h) = —-—, h > 0, (30) 


with 


Q = max (/ijvs, 75 - (ps + ws)). 
Using the fact that 'Jg — f-i^G = — 1), we dehne 


(31) 


4>G{h) = < 


gMG(^G-l)^ _ X 


Hg > 0 . 


eich 

-, /iG = 0, h > 0. 


(32) 


V 7 G 

This scheme is positively invariant and is qualitatively stable, which means that it has the 
same equilibria than system ([I])-(I2]), and the stability/instability properties of the equilibria 
are preserved, at least locally, whatever the stepsize h > 0 , [52] . 


4.2 Parameters ranges and ecological zones of biomass produc¬ 
tions 

To provide relevant numerical simulations, one need to use ecologically meaningful param¬ 
eters ranges and values. Thus, after extensive literature review, we found the following 
parameters ranges: 
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Table 2: Parameters values found in literature 


Parameters 

values 

References 

/ (1/^) 

0 - 1 

Langevelde et al. (2003) |3T] 


0-2 

Accatino et al. (2010) 

7g 

0.4(1) _ 4_g(2) 

(1) Penning de Vries (1982) [58] 

(^) Menaut et al. (1979) 

7s + 7Afs 

0.456 - 7.2 

Breman et al. (1995) [15] 

yNs 

0.03 - 0.3 

Accatino et al. (2010) [3] 


0.4 

Langevelde et al. (2003) [3T] 

hs 

0 - 0.3 

Langevelde et al. (2003) [8T] 

he 

0 - 0.6 

Langevelde et al. (2003) [8T] 

Vg 

o.i(“) -1(^) 

(“) Van de Vijver (1999) |1H] 

(^) Accatino et al. (2010) [5] 


0.2 - 1 

Abbadie et al. (2006) [T] 

Vs 

0.02 - 0.6 

Accatino et al. (2010) [3] 


0.66 

Reinterpretation of Gignoux (1994) [28]. 
Reinterpretation of Langevelde et al. (2003) |3T] 

Us 

0.05 - 0.2 

Wakeling et al. (2011) [H] 


As stated previously, we will focus our numerical simulations on three ecological zones 
of the African continent 

Region 1 is a semiarid area with mean annual rainfall comprised between 300 mm.yr~^ and 650 
mm.yr~^ where there is a low biomass production and few hre occurrence (says one 
hre event every ten years) if any. 

Region 2 is a mesic area with a mean annual rainfall that varies between 650 mm.yr~^ and 1100 
mm.yr~^ which is an intermediate biomass production zone and where we can have 
on average one hre event every four or hve years and sometimes less. 

Region 3 is a high biomass production zone, in which we can have one or two hre events per 
year, i.e. a humid tropical area with a mean annual rainfall between 1100 mm.yr~^ 
and 1800 mm.yr~^. 

Our aim is to assess the diherent outcomes of the IFAC model along with the inhuence of the 
variations of ac, ctns (for which there is no direct information in the published literature) 
and the hre period r. The parameter ranges in each of the regions are summarized in Table 

El 
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Table 3: Parameters ranges in the three ecological regions 


Parameter 

Region 1 

Region 2 

Region 3 

r {yr) 

> 5 

2-5 

0.5 - 3 

Kt {t.ha~^) 

30 

80 - 90 

110 - 120 

Kg {t-ha~^) 

0-5 

5-10 

10 - 20 

7g {yr~^) 

0.4 - 2 

2 - 3.5 

3.5 - 4.6 

75 {yr~^) 

0.2 - 0.8 

0.2 - 1 

1.5 - 2.7 

7NS {yr~^) 

0.256 - 1.2 

1.2 - 2.5 

2.5 - 4.5 

yNs iyr~^) 

0.1 - 0.25 

0.07-0.1 

0.02 - 0.07 


In addition, reinterpreting experiments that concern Region 1 and Region 2 and reported 
in February et ah (2013) [2S] we derived ac (in ha.t~^.yr~^): 0.1843 - 0.9984 for Region 
1 and ac'- 0.2470 ~ 1.6287 for Region 2. Moreover, several studies located under different 
rainfall compared grass production under and outside a tree crown. A synthesis was proposed 
by Mordelet & Le Roux (see Abadie et al. (2006) Page 156 [1]) that emphasized that the 
relative production (within to outside) is a decreasing function along the rainfall gradient. 
We re-interpreted the results as to derive reasonable values for ai\fs in each of the three 
regions, using the subsequent reasoning. 

Assuming that the measurements were made, free of hres, in grass stands having reached 
equilibrium, and letting Gu and Go be the equilibrium values under and outside crown, 
respectively. We can write according to the model: 


G„ = 


(33) 


The ratio considered by Mordelet & Le Roux (1995) |36] 
156 m) is: 


(see also Abadie et al. (2006) Page 


(34) 


i.e. the ratio of grass production under and outside a tree crown. Assuming yc 
simplify as: 


(5g = 1 - 


o'nsT 

IG 


(^NS 


(1 - 5g)ig 

f 


0 we can 
(35) 


T is the woody biomass density to be computed at the scale of an isolated, full grown tree 
(having reached the maximal height considering the local climate) in any of the three regions. 
We propose to relate T to Kt as 

~ X Kt 
S ^ 

where S is the woody cover characterizing the maximal density Kt and e g] 0, 1[ is a coeffi¬ 
cient expressing that an isolated tree has less influence on grass production that a complete, 
closed canopy stand corresponding to Kt- 
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We used the value of Sq = 1-58 (resp. 1.25 and 0.75, 0.75) of Mordelet & Le Roux (1995) 
[5B] (see also Abadie et al. (2006) Page 156 [T]) that corresponds to Region 1 (semiarid 
region) (resp. Region 2, Region 3). Note that 6g values above 1 express a facilitative 
effect of trees for grass while values below correspond to a depressing effect. Using also the 
estimated values for Kt and 7 g, we deduce the ranges of variation of a^s (in ha.t~^.yr~^), 
summarized in Table 01 

Table 4: Variation range of aisfs in Region 1, Region 2 and Region 3 following re¬ 
interpretation of Mordelet & Le Roux results 




7g 

Kt 

£ (min) 

£ (max) 

5 

airs (min) 

aNS (max) 

Region 1 

1.58 

0.6 

30 

0.4 

0.75 

1 

-0.029 

-0.0155 

Region 2 

1.25 

2.8 

85 

0.2 

0.67 

1 

-0.0412 

-0.0123 

0.75 

0.0123 

0.0412 

Region 3 

0.75 

4.2 

115 

0.1 

0.15 

1 

0.0609 

0.0913 


Note that the range of values for a^s is due to a large uncertainty on e. Moreover as a 
straightforward consequence of the results of the reference study (Modelet & Le Roux (1995) 
[36] . Abadie et al. (2006)) a^s is likely to be negative in Region 1 (shading improves grass 
production in Region 1, which is also in good agreement with results of Moustakas et al. 
(2013) [37], Belsky et al. (1989) [13], Weltzin & Coughenour (1990) [57]). 


5 Numerical simulations and discussion 

5.1 Results for Region 1 

In semiarid areas, the main mechanisms that govern the ecological processes include 

(Ml) water limitation on tree growth (Baudena et al. (2014) [TT]) and on grass biomass 
standing crop (Kg) 

(M2) tree - grass competition, which has an especially strong competitive impact on tree 
seedlings (February et al. (2013) [25] I 

(M3) unfrequent hre may reduces woody cover and grass cover but grass biomass recover 
quickly after hre (Baudena et al. (2014) [TT]) while low values of grass biomass limit 
the impact of competition on sensitive tree biomass. 

Point (M2) suggests that the sensitive tree vs. grass competition parameter ug has relatively 
large values while point (Ml) along with point (M3) suggest that woody cover is controlled 
principally by water availability and secondarily by unfrequent hres (Sankaran et al. (2005) 
[10]). Thus the depressive effect of scattered woody cover on grass biomass mainly results 
from reduced light availability and root competition for soil water (see Walker et al. (1981) 
[50]L Nevertheless, somme references (Belsky et al. (1989) [13], Weltzin & Coughenour 
(1990) [51]) also emphasized the facilitation role of scattered tree on grass biomass in East 
African semiarid savannas. Indeed, compared with the open situation, the highest grass 
production was recorded under acacia and baobab trees (Mordelet & Menaut (1995) [35] . 
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Belsky et al. (1989) [13], Weltzin & Coughenour (1990) [5l]), which are known to have a 
low light interception and only induce a slight limitation to photosynthesis while shading 
improves the water balance under the canopy (Barbier et al. 2008) |H]. This finding of 
Mordelet & Menaut (1995) [36|, Belsky et ah (1989) [I3], Weltzin & Coughenour (1990) 
[5T| can also be explained by the soil enrichment by nitrogen fixing species, like acacias trees 
which results in a yield increase. Therefore, in semiarid areas the main ecological vegetation 
types that are observed depending on annual rainfall (which only varies tree/grass ratio) are 
savannas (February et al. (2013) [2S], Baudena et al. (2014) [TT], Accatino et al. (2010) [5]) 
and sometimes forest in the sense of low dry forest, and/or thickets (Walker et al. (1981) 
[50] , Couteron & Kokou (1997) [IZ])- The outcome of the IFAC model in semiarid areas is 
also in adequacy with this previous features. Indeed, let us consider the following table of 
parameters values 


Table 5: Array of parameters’ values for Region 1 


75 INS IG 

^^NS 

us 

hs 

Vs 

0.3 1 0.6 

0.15 

0.1 

0.2 

0.5 

II 

-<i 

II 

CO 

II 


For values in Table lO] we compute 



nl 

3.2222 

2 


and we derive figure [H 
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(B) Level curve of the threshold in Region 1 



Figure 2: Level curve of the threshold illustrating that as rjc, hg and aNs increase, 
decreases and system ([ID-([2]) is liable to move from a savanna/grassland state to a forest 
state or to a multistability involving the forest solution. Recall that the forest solution is 
stable (resp. unstable) whenever p^ is lower (resp. greater) than unity. In (A) pc = 0.3, in 
(B) pg=0.5. 


Recall that according to relation (|9]), 

0 , f 0.8775 for pc = 0.3 

^ ^ hG < I g 5034 ^ 0 5 

Moreover, together with our data estimation of ac and utvs in Region 1, we found that 
77^ > 1. Setting pc = 0.3, rjc = 0.6, and 

• when ac G [0.92, 0.95[, we have 77^ > 1 and Pt < 1 

• when (Tg e [0.95, 0.9984[, we have 77^ < 1 and pr < 1. 

Finally, when pc = 0.5 one gets pt > 1- Therefore, hgure [2] together with the previous 
discussion on TZq and px illustrates either case 1 to case 6 of Table 13.41 Thus, when external 
disturbances (such as herbivory) on grass biomass are low and ug has relatively large values, 
the IFAC model predicts either a stable savanna state, a stable forest state, a stable grassland 
state or a multistability involving savanna and/or forest and/or grassland (see also case 1, 
case 2, case 3, case 5 and case 6 of Table iRTjl . Furthermore, when external disturbances 
on grass biomass become more important, the grassland solution becomes unstable and the 
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IFAC model predicts either a stable savanna state, a stable forest state or a bistability 
involving savanna and foret states (see also case 1 and case 4 of Table 123]) • Consequently, 
one can observe that the non sensitive tree vs. grass interaction parameter a^s and the 
additional death rate of grass biomass due to external disturbances pc are likely to be 
influential on the IFAC outcomes in Region 1. Nevertheless, with a mean annual rainfall of 
300-400 mm there is a wide array of references evidencing the probable bistability of desert 
(bare soil) and thickets (in the African Sahel) (see Couteron & Kokou (1997) [IT], Lefever 
et ah (2009) [32], Barbier et ah (2008) [9]) or desert and grass (Namibia) (see Tschinkel 
(2012) [HI, Fernandez-Oto et ah (2014) [23]b 

5.2 Results for Region 2 

In Region 2 which corresponds to a mesic area, the main mechanisms that regulate tree-grass 
interactions also include mechanisms (Ml) and (M2) stated previously for semiarid areas. 
In addition to (Ml) and (M2), in mesic areas, hres are more frequent than in semiarid areas 
since water availability favor grass biomass production which constitutes the fuel for hres 
(we denote this new mechanism (M4)). Grass-hre feedback (mechanism M4) together with 
mechanisms (Ml) and (M2), maintain both forest and savanna occurrences in mesic areas 
(Baudena et al. (2014) [TT]'). Indeed, grasses beneht from the openness of the landscape 
after hres, since they recover faster than trees seedlings, thus determining a positive feedback 
mechanism that enhances savanna presence. The IFAC model also predict a shift from a 
forest state to a savanna state as a^s decreases and/or when ac increases, which also agree 
with Baudena et al. (2014) [TT] results. Therefore, savanna and forest vegetation types 
clearly appear as alternatively stable states as found by Staver et al. (2011) [33], Staver and 
Levin (2012) ]33], in our case depending on ac and a^s variation in Region 2. 

Let us consider the following table of parameters values 


Table 6: Array of parameters’ values for Region 2 


7S INS IG 


ojs 

ris 

VG 

hs 

0.4 2 2.8 

0.08 

0.1 

0.5 

0.6 

0.1 

Kt=85, Kg = 7, pc = 0.3 


Using parameters values in Table [O] one has: 




14.5 

14 


and we also derive hgure [3] and hgure 01 
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Figure 3: Level curve of the threshold illustrating that system ([I])-® is liable to move 
from a savanna/grassland state to a forest state or to a multistability involving the forest 
solution together with r and a^s variations. Recall that the forest solution is stable (resp. 
unstable) whenever is lower (resp. greater) than unity. 


(A) Level curve of the threshold p.^. in Region 2 



(B) Level curve of the threshold p in Region 2 



Figure 4: Level curve of the threshold px illustrating that system (II])-(I2]) is liable to move 
from a savanna/forest state to a grassland state or to a multistability involving the grassland 
solution in relation to r and ac variations. Recall that the grassland solution is stable (resp. 
unstable) whenever px is lower (resp. greater) than unity. In (A) pc = 0.2, in (B) pc = 0.3. 
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Since for parameters values in Table [6] one has > 1 then, 


Pg>1 


T > 


0.3524 for /xg = 0.2 
0.3665 for/rG = 0.3. 


Moreover, for our estimation of aNS and ac one also has > 1 and TZq > 1. Therefore, 
hgurelHand figure [3] illustrate, either case 1, case 2, case 4 or case 5 of Table ITTI 

In summary, the parameters that are likely to be influential on the IFAC outcomes in 
Region 2 are the external disturbances on grass biomass parameter /ig, the grass vs. sensitive 
tree parameter ctg, the non sensitive tree vs. grass interaction parameter a ns- In addition 
to that previous parameters, one can also mention the hre return time r. In other words, 
the previous analysis reveals that in Region 2, in addition of stability of forest, stability of 
savanna and bistability of forest and savanna as in Region 1, one can observe stability of 
grassland and also multistabilty situations involving grassland solution with relatively low 
values of ac in comparison with Region 1. 


5.3 Results for Region 3 

Region 3 corresponds to humid tropical areas where rainfall availability favors biomass pro¬ 
duction of both woody and grasses components. The grass-hre feedback possibly leads to 
stability of either savanna or forest in Region 3 depending on hre return time. Indeed, grass, 
particularly abundant in these wet areas, becomes an extremely good fuel in the dry sea¬ 
son, which promotes hre occurrence and increases hre intensity and impact (Baudena et ah 
(2014) dH, Higgins et al. (2008) [29]). When the hre return time is large, the trees have 
the time to grow above the hame zone and to reach canopy closure and then outcompete 
grasses. Therefore, relatively large return time favor forest state in Region 3 (Staver and 
Levin (2012) [15 )• Moreover, if the hre return time is small then trees don’t have the time to 
reach canopy closure and therefore let grasses which regrow quickly in the open space after 
hres form either a stable savanna state or a stable grassland state. We illustrate hereafter 
that these features are also predicted by the IFAC model. 

Consider 


Table 7: Array of parameters’ values for Region 3 


IS 

INS 7G ^J-NS 

Us 

75 

7G 

h5 

FG 

2 

3 4.2 0.06 

0.1 

0.5 

0.6 

0.1 

0.2 

Ar=115, Kq = 15 


In this section, we will refer to a particular area, namely the Lamto region in Ivory Coast 
(see Menaut et al. (1979) [35], Mordelet & Menaut (1995) [36] )• Thanks to Abadie et al. 
(2006) Page 156 [T] (see also Mordelet & Menaut (1995) [36]) data report of grass biomass in 
the canopy and open situations in Lamto and by reinterpreting their results we have derived 
the range 0.0609 < a^s < 0.0913. 

Using parameters values in Table H one has: 



nl 

35 
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We also derived the following figure |5l 



Figure 5: Level curve of the threshold px illustrating that system ([I])-® is liable to move 
from a savanna/forest state to a grassland state or to a multistability involving the grassland 
solution in relation to r and ac variations. Recall that the grassland solution is stable (resp. 
unstable) whenever pT is lower (resp. greater) than unity. 

Since for parameters values in Table [7] and according to our estimation of a^s one has 

• TZ^ < 1 , < 1 , 

• > 1 then, 

> 1 r > 0.2291. 

Values of r are not expected to be under 0.5 (i.e. 2 hres per year), a minimum which 
corresponds to sub-equatorial climates with two dry seasons. We therefore consider 
that this condition is always fulhlled. 

Therefore, one can deduce that hgure |5] illustrates either case 8, 9 or case 10 of Table iTH 

The previous analysis highlighted the importance of the grass vs. sensitive tree compe¬ 
tition parameter, (Tg, the non sensitive tree vs. grass competition parameter, cttvs, and the 
hre return time, r, in controlling the outcomes of the IFAC model. Comparing to results of 
Region 2 (in terms of having the forest, the grassland or the savanna as reachable solution), 
one note that the IFAC model fairly has the same outcomes in Region 3 as in Region 2. 
Figure [5] further show that for low values of ug, say less than 0.6, r seems to have very 
limited influence. It is stronger for large values for which increases in r make pt decreases 
under 1 and therefore destabilize the grassland solution. 

In summary, in Region 3, we observed that the hre return time along with the grass vs. 
sensitive tree competition parameter ac and the non sensitive tree vs. grass competition 
parameter a^s strongly inhuence the outcome of the IFAC model. Indeed, depending on 
these parameters variations and values, the IFAC can converge either to a grassland state, 
to a savanna state, to a forest state or to a multistability involving forest state and either 
grassland or savanna while environmental conditions in this Region would systematically 
allow forests in the absence of hre (Staver and Levin (2012) |11]). Therefore to favor a forest 
state in Region 3, one could implement policies in order to have relatively large hre return 
time (says a hre frequency, f, lower than one hre per year: / < 1) and vice versa if tracts of 
savanna are to be kept against forest encroachment as habitats of large grazing mammals. 
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6 Conclusion 


In this work, we presented and analyzed a new kind of mathematical model for tree-grass 
interactions in savanna ecosystems either fire prone or not. It is an extension of a continuous¬ 
time model, called the COFAC model, studied in Yatat et ah (2014) |52]. The model 
presented here, that we call the IFAC model, is based on a system featuring impulsive 
differential equations and thereby aims to acknowledge the discrete nature of fire events. 
The analytical study of the IFAC reveals a desert equilibrium, a forest equilibrium and two 
periodic solutions: the grassland periodic solution and the savanna periodic solution. The 
analytical study also reveals seven ecological thresholds (7^^, 71%, ^ti Pti Pt)- 

These thresholds dehne in parameter space regions of monostability, bistability also found 
with the models of Accatino et al. (2010) [3], De Michele et al. (2011) [T3], Yatat et ah 
(2014) 152] ■ They also dehne regions of tristability as in Yatat et al. (2014) [52] with respect 
to the equilibria (desert and forest) and periodic solutions (grassland and savanna). The 
specihcity of IFAC is to also present periodic behaviors which depict huctuations in woody 
and/or grassy biomass and cannot be yielded by fully continuous time models. Therefore, in 
case of transition from a vegetation type (forest, savanna or grassland) to another, the change 
is done progressively which is more ecologically meaningful than the abrupt changes observed 
for hres-continuous models such as in Accatino et al. (2010) |3], De Michele et al. (2011) 
[19] and Yatat et al. (2014) [52]. Such abrupt changes have been criticized as unrealistic 
by Accatino et al. (2013) [2] (see also Beckage et al. (2011) [12]) who advocated stochastic 
models in lieu of continuous time formulation of hre impact on vegetation. The present 
approach however demonstrates that a more realistic modeling of hre can be introduced 
within the framework of continuous time models while keeping the potential for analytical 
exploration of the main outcomes of the model. Something which is not possible with the 
aforementioned fully stochastic models. 

As in Yatat et al. (2014) [52], we found that the competition parameters (Tq which ex¬ 
presses the asymmetric competition exerted by grasses on sensitive trees (shading and soil 
resource preemption) and ajsis that expresses the asymmetric competition of non sensitive 
trees on grasses (shading and soil resource preemption) are bifurcation parameters of the 
IFAC model along with hre frequency (which strongly inhuence the convergence outcomes 
of the IFAC model). The analytical study of the IFAC model also reveals three particular 
values r*, and aQ (see relation fl53]l . Page|12]) that delimit regions of stability/instability 
or forest and grassland solutions in relation to r, ac and ap^s respectively (see relations f[T3]) . 
Page El and flTT]) . Page|9]). Moreover, considering three ecological biomass production zones 
indexed by hres frequency and by carrying capacities of both trees and grasses biomass, 
allows us to point out several scenarios for IFAC convergence that depend on ctg, cns and 
r values and variations. These outcomes of the IFAC model are qualitatively in agreement 
with results of Baudena et al. (2014) [TT], Staver et al. (2011) [13], February et al. (2013) 
[25] . Accatino et al. (2010) [3], Couteron & Kokou (1997) [IT], Staver & Levin (2012) [41] . 
Distinguishing three ecological zones allowed us to verify in which contexts the possible bifur¬ 
cation parameters are actually inhuential or not. This analysis highlighted the pervasiveness 
of crv 5 (he. the depressive or facilitation effect of grown-up trees on grasses) in all the three 
zones. It also emphasized the influence of ac (i-e. depressive effect of grasses on small trees) 
in the two zones (2 and 3) with sufficient rainfall to allow medium to high grass production. 
In these two zones, and especially in zone 3, the fire return time (r) appeared also influential. 
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As already mentioned by Yatat et al. (2014) [52], and verified here for the IFAC model, the 
competition parameters ac and a^s which embody direct tree-grass interactions deserve 
an increased interest and should be the focus of adhoc observations and experiments as to 
better assess their ranges of variation in the different ecological regions. 

Although the IFAC model presented in this work and the COFAC model presented in 
Yatat et al. (2014) [52| qualitatively display strong similarities, the IFAC model is richer. 
Indeed modelling fire events as pulse phenomena leads to a relaxation of stability conditions 
of both forest and grassland solutions and it increases parameters ranges for which bistability 
situations involving forest, grassland and savanna can occur. This particular property of the 
IFAC model may explain, along with bifurcation parameters and its periodic outcomes, 
many changes in tree-grass interactions in fire-prone ecosystems. Thanks to this particular 
property, the IFAC model, which moreover displays periodic outcomes and bifurcations 
according to well-identified parameters, is able to account for many dynamical scenarios 
observed in savanna-like ecosystems from the fringes of the desert to the boundary of the 
wet forest. 
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Appendix A: Proof of Lemma [3 

Defining 

Ts{t) = x{t), 

TNsit) = y{t), 

G{t) = G*{t)+z{t), 


(36) 


where x{t),y{t) and z{t) are small perturbations. Every solution of the linearized equations 
can be written as 




a;(0) 


I y{t) I = m I 1/(0) 


\ 

Here $ is a fundamental matrix and satisfies, 

d<h(f) 


(37) 


dt 


= DF{0;0;G*{t))m) 

^ 7s — (hs + + crGG*(f)) 

V « 


Ins 

—y-Ns 


0 

0 




7g 


-aNsG*{t) '^G —— hg 

Ag 


m 




m 

and *F(0) = Id^s. Moreover the resetting impulsive condition of system (II])-(I2]) becomes, 

(39) 


x(nr+) 

//(nr+) 

z{nT^) 


1 — r/5ta(G*(r)) 0 0 

0 1 0 
0 0 1 — 7/g 


x{nT) 

yinr) 

zinr) 


A monodromy matrix M of system (ITl) — (j2 |) , is: 

1 — r/5tc(G*(r)) 0 0 

M= I 0 10 

0 0 1 — r/G 


<l>(r). 


with 


= exp DF{0;0;G*{s))ds 


(40) 


(41) 


Moreover using Lemma [21 a direct computations leads 

0 

DF{0-,0-,G*{s))ds = I DF® 0 

0 DFG) DF^^'^ 
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where 

L)F(i) = ( 7 s - (/^s + -(Tg [ G*{t)dt, 

= 'fNST, 

DF(3) = UsT, 

DFG) = ^ 

= -aNS j G*{t)dt, 

DF^^) = -;xG(7^^-l)r-21n(l-77G)• 

Consider the sub-matrix B dehned as follow: 

/ DFd) L)F( 2 ) \ 

\ ) ■ 

Recall that eigenvalues of the matrix B are root of the quadratic equation 

— frace(B)A -|- det(B) = 0 


( 42 ) 


(43) 


and to characterize real part of eigenvalues of matrix B, following Routh-Hurwitz cri- 
terium (see Section 1.3.5 Page 72 of Augier et al. (2010) [7]), one need only to study the 
sign of frace(B) and det(B). Let 

A = tr{B) = 'ysr ’ (44) 


where 


n 


Thus, if 77 < 1 then Al < 0. 
Moreover, let 


_ Js _ 

l^s + + I^NS + O-cGint 


> 0 . 


^ det(B) ^ rfij^s {(ps + + (JcGint) {l - TIq) > 


(45) 


where 



_ 'yS^^NS + _ 

fJ^Nsid'S + ^s) + (^G^-NsGint 


> 0 . 


Thus, if 77^ < 1 then S > 0. 

Moreover, one also has 77 < nl. 

Therefore, if 77^ < 1 then s(B) < 0, where s denotes the stability modulus (i.e. the 
maximum of the real part of eigenvalues). 

From expressions flTn|l and flTT]l we deduced that eigenvalues (^i, ^2 and .^3 of the monodromy 
matrix M are 


6 

6 


(1 - r]sw{G*{T)))e^\ 

1 - hG 


( 46 ) 


where Ai, A 2 G sp(B). 

Since 0 < 1 — rjsw{G*{T)) < 1, if 77^ < 1, then 0 < .^1 < 1 and 0 < .^2 < 1- 
Moreover, since > 1 then ^3 < 1. Indeed, 


^3 < 1 < 1 - pg 77 1 < (1 - ^ > 1 . 
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Finally we deduce that the grassland periodic solution Eq = (0;0;G'*(t)) is locally asymp¬ 
totically stable if < 1 or (TZq > 1 and < 1), is locally stable if (TZq > 1 and pt = 1) 
and is unstable if (TZq > 1 and px > !)• This ends the proof. 


Appendix B: Proof of Theorem [T] 


* Case 1: pc > 0. 

Solution G of system ([I])-© satisfy 

G'(i) < 70 (l - G(«) 

= (l-hG)G(4)- 

From Lemma 1.3 page 15 in [8] we deduce that 


(47) 


G(i) < G(0) ( n (70 fl - 


K 


Thus, for TZq < 1 we have lim G(t) = 0 and solutions Ts and T^s of system ([I])-([2]) satisfy 

i—4-I-00 


T's = {'ysTs + InsTns) ( 1 - ^ k ^ ~ Ts{ps + ^s), 


'^Ns — ^sTs — PnsTns- 


(48) 


System flT8)l does not admit periodic solution (see Appendix B in Yatat et al. (2014) [52]), 
thus using the Jacobian matrix of system (ITS]) we deduce that 

• if < 1 then, {Ts,Tns) -t (0,0), 

• if TZx > 1 then, (Ts,Tns) —t (Ts,Tpfs)i where (Ts,Tpfs) given in (nai. 

At the end, we deduce that if 7Z^ < 1 and TZ^q < 1 then, the desert equilibrium Eq is GAS 
i.e., point 1 of Theorem [T] holds. The forest equilibrium Ex is GAS whenever TZ^ > 1 and 
TZq < 1 i.e., point 2 of Theorem [1] holds. 

Now suppose that TZx < 1 and TZq > 1. Solutions Ts and T^s of system ([I])-(j2|) satisfy 


dt 

dTp^s 

dt 


< ilsTs + 7nsTns) 1 


< ^sTs — t<^NsTNS, 


Ts + Ti 


NS 




Tsips + ^s), 


t^tk 


Ts{tl) < Ts{tk), 

TNsit'l^) < TNsitk), t = tk 4+1 = tk + T. 
Let consider the upper system 


(49) 


(50) 


^ = hsU + 'JNSV) ( 1 


dv 

— = UJsU — Pns'v, 


U + V 

Kx 


u{ps + UJs), 


(51) 
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Since TZ^ < 1, {u{t),v{t)) —?• (0,0). Thus {Ts{t),TNs{t)) (0,0). Furthermore, solution G 
of system ([I])-([2]) admits as limiting system 

= 7G(l-£)G-f,cG (#4. , 52 , 

I ^i^k) ~ (1 ~ hG)G^(4) t = tk- 

System fl52|l admits at most two solutions: the trivial solution, 0, which always exists and 
the periodic solution G*{t) which is ecologically meaningful if > 1 where G*{t) is given 
by ([8]). Now we turn to check stability results of solutions of system fl52|l through small 
perturbations approach and Floquet’s theory. 

• Setting G{t) = x{t) where a; is a small perturbation and verify x{t) = (j){t)xo, where 0 
verify 

0'(^) = hG (^G - 1) 

and 0(0) = 1. The resulting impulsive condition becomes 

x{nT~^) = (1 — riG)x{nT). 

Following the Floquet’s theory, the zero equilibrium is locally asymptotically stable if 

Ao = (l-pG)e''«('"«-')"<l. 

Since 

TIq > 1 , Aq < 1 if and only if (1 — < 1, i.e. < 1. 

Moreover, for < 1 the positive solution G*{t) is undehned then the desert equi¬ 
librium is globally asymptotically stable. Finally, we deduce that the desert solution 
(0, 0, 0) is globally asymptotically stable whenever 71^ < 1, TZq > 1 and p^ < 1. Point 
3 of Theorem [T] holds. 

• Now, setting G{t) = G*{t) + x{t) where x is a small perturbation and verify x{t) = 
0 (f)xo, where 0 verify 




mg {K - 1) - ^G-(i) 


0(f) 


and 0(0) = 1. The resulting impulsive condition becomes 

x{nT'^) = (1 — T]G)x{nT). 

According to the Floquet’s theory, solution G*{t) is locally asymptotically stable if 


Ag* = (1 - Vg) exp < pg ("^g - l) t - 


27g 

Kg 


“(n+l)r 


G*{t)dt ^ < 1. 


Following Lemma |2], 

f{n+l)T 


G*{s)ds = — {ln(l-P g)+ hG (^G - l) 
7g 


(63) 
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Thus 


(54) 


Ag* = (1-r7G)exp{-/XG(^G - l)^-21n(l-r^c)} 

= exp{-/iG(7^G--77 g)}- 

Since > 1, we have: 

Pg > 1 ^ (1 - Pg) exp{pG(^G - 1)^} > 1 

ln(l - t]g) > -Pg(^g - 1)'^ (55) 

-Pg(^g - 1)^ - - Pg) < 0 


and we deduce 


Ag* ^ 1- 


Thus solution G*{t) of fl5^ is globally asymptotically stable because the zero solu¬ 
tion, in this case, is unstable. Finally, we deduce that the grassland periodic solution 
(0, 0, (j'*(t)) is globally asymptotically stable whenever 71^ < 1, TZq > 1 and > 1. 
Point 4 of Theorem [T] holds. 


-k Case 2: pc = 0. 

The proof of points (i) and (ii) of Theorem [T] is fairly the same as the proof of points 3 and 
4. Indeed we first set, only in System fl52|l . /ig = 0 and next, we substitute pciJ^G ~ 1) ^y 
7 g in the rest of the proof. 

Solution G of system ([I])-([2]) satisfy 


dG 

dt 

Gitt\ 


= (1 - r]G)G{tk) t = tk. 


(56) 


Since = (1 — pg) exp{ 7 Gr} < 1, it follows that G(t) —> 0. Therefore, solutions Ts and 
Tns of dll)-)© satisfy system (145]) . Since TZ^ < 1, one has {Ts, T^s) —> {Ts, Tvs)- Point 
(hi) of Theorem [T] holds. 


Appendix C: Proof of Theorem [2] 

Taking new variables Ts{t) = Tjvs(t) = G{t) = then system ([T]) 

becomes. 


x{t) 

= -UJS - Ps- o-gc^ + (7s + 7ivse^e 

yit) 

= -PNS + OJsG'e-y, tn+l = tn + T, 
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) f 7^ fn- 


(57) 


Let X = C^([0, r], M^), Z = C^([0, r], M^) x C^([0, r], M^) and for u = {x, y, z) G X, 


lull = max |a;(f)| -|- max |p(f)| -|- max \z{t)\. 

te[o,r] te[o,r] te[o,r] 
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Then X, Z are Banach spaces when they are endowed with the above norm 
Let, 


L : Dom{L) G X ^ Z, 


and 


where 




/ Ax{tn) 
Ay{t. 

\ Az{t 




-Us - ys- + (7s + 7 vse^e 1 - 
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A direct computation leads to 


KerL = 


1 - — ) - aNse^ - iiG 
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and 


ImL = 




G Z 


t e [0, r] 


fj l{t)dt + a = 0 
JJ" m{t)dt + 6 = 0 
n{t)dt + c = 0 


Since ImL is closed in Z, L is a Fredholm mapping of index zero. Indeed, 

Index{L) = dim{KerL)—dim{CoKerL) = dim{KerL) — {dim{Z)—dim{ImL)) = 3—(6—3) = 0. 

Thus following (Gaines and Mawhin (1977) [27], Page 12), there exist two continuous pro- 

P ^ L_ „ Q 


jectors P and Q such that the sequel X —)■ DomL 
and KerQ = ImL = Im{I — Q)- It suffices to choose 


and Q 




One can verify that LP 



= Ox and QL 



Z is exact i.e ImP = KerL 


l{s)dt + a 
JJ" m{s)dt + b 
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Furthermore, the generalized inverse Kp : ImL —)• KerPnDom{L) of the map L : KerPn 
Dom{L) —>■ ImL is given by 

Jq Z(s)ds + a 


Kr 




= I Jq m{s)ds + 6 
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Indeed, let u = e KerP fl Dom{L), {g,r) = {{ 91 , 92 , 93 ), {ri,r 2 ,r 3 )) e ImL, we 

have 

KpL{u{t)) = Kp{u, Am) 

= 'u(s)(is + Am 

= u{t) — m(0) + m(0) — u{t) (58) 

= u{t) — P{u) 

= u{t), because u G KerP 

LKp{g{t),r) = L{jl g{s)ds + r) 

= {9{t), - lo 9{t)dt) (59) 

= {g{t), r) because {g,r) G ImL. 


Thus, 
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(63) 

Clearly, QN and Kp{I — Q)N are continuous then for any open bounded set f2 C X, 
QN{pL) is bounded. Furthermore, let ti, t 2 G [0, r], u{t) = {x,y^z){t), 
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rti 


f{s,u{s))ds- / f{s,u{s))ds 


T 


f{s, u{s))ds + a 


T 

ri2 


/(s, u{s))ds + a 
(^2 ~ ti) 


f{s,u{s))ds 


f{s, u{s))ds + a 


\t2 — tl 


< \t 2 - til max |/(t,M(t))| + 

te[o, r] r 


r max |/(t, u{t)) \ + a 
te[o, r] 


< |t 2 - til 2 max |/(t,M(t))| +- 

iG[ 0 , t] T 


(64) 


and 


\Kp{I — Q)N{u{t))\ < |a| + r max |/(t,M(t))| + |a| + r max |/(t,M(t))| 

te[o, t] tG[ 0 , t] 

< 2 ( |a| + r max |/(t,n(t))| ) . 

V i6[0, d J 


( 66 ) 


Then using relations fIMl) . fl65|) and the Arzela-Ascoli’s theorem (Sonntag (1997) [12] Theo¬ 
rem 3.1, Page 314) we deduce that 
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Kp{I — Q)N{Q) is compact. Thus, is a L—compact mapping on The isomorphism J 
of ImQ onto KerL may be dehned by 


J : ImQ —>■ X, 





Now we reach the position to search for an appropriate open, bounded subset for the 
application of the continuation theorem, i.e we search Mq such that every r-periodic solution 
of system ([I]) — ([2]) satished |a;(f)| + \y{t) \ + \z{t)\ < Mq with 0 < t < r. 

Corresponding to the operator equation Lx = (3Nx, (3 G (0,1), we have 


x{t) 

y{t) 

z{t) 


= P 
= 

= /S 


(^s - Ls - o-Ge"" + (7s + 'iNse^e-Q ( 1 
^TVS - 

7g ( 1 


K. 


t ^ tn- 


aNse^ - Lg 


= 13 [-fiNS + ujse^e , t^+i = tn + T, 

x{t~^) — x{t) = /3\ia{l — r]sw{eQ), t = t 
y{t^)-y{t) = 0 , n = 0 ,l, 2 ,..., 

^ z{t^) -z{t) = /?ln(l-7G)- 


Suppose that {x{t),y{t), z{t)) G X is an arbitrary solution of system 
(3 G (0,1). Integrating on both sides of fl^ over the interval [0, r], we obtain 


( 66 ) 

for a certain 


Is 


K, 


{e^ + e^) + 'yNse^e~ 


1 - 


Kt 


- (JGe 


dt = {us + /is - 7 s)t 

— ln(l — ?75t(;(e^('^))), 
use^e~^dt = Uns'^, 


7g 

0 iKg 


+ CTvse^ 


dt = (7g -/ic)^ + ln(l - 7 g)- 

(67) 


Note that assumptions of Theorem [2] lead 

(7g - Lg)^ + ln(l - Vg) > 0. 

Since X is a Banach space and {x{t),y{t), z{t)) G X, there exist fj, r/, r and r such 

that 


xil) = max xit), xil) = min xit), 

0<t<T - 0<t<T 

y{v) = max y{t), y{r]) = min y{t), 

0<t<T — U<I<T 

z{t) = max z{t), z{t) = min z{t). 

0<t<T 0<t<T 


( 68 ) 


It follows from system fl67)) that 


VI 

i^s 

+ ys)T 


-(TgC^ 

+ 

{is + iNse^e Q A 

< 

i^s 

+ ys)r 

+ acTe^^'^'’ + 

L 

T ^ 

(7s + lNse^e~Qdt 

< 

i^s 

+ ys)r 

+ + 

1 

T 

{is + -iNse^)dt 

< 

i^s 

+ /is + 

7s)^ d 

acTe^ 


+ lNs'^ey^^\ 
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_j_ 

Kt 


dt 


(69) 




















and 


Recall that 


Since 


then 


Moreover, from 


we deduce 


/ \y{t)\dt < IINS'T+ / ^\dt 

'o Jo 

< 2fi]S[sT 


\z{t)\dt < (7G + hG)^ + y 

< ilG + yG)r + (7g - yG)r + ln(l - ^g) 

< 27 Gr + ln(l - 7 g)- 


Pg > 1 (7g - PG)'r + ln(l - 7 g) > 0 

27 Gr + ln(l - rjc) > 0. 


< (7g - Pg) 


r — / z(t)dt 

Jo 

< (7g - Pg)^ + ln(l - 7 g), 


,, ^ , ln(l -7 g) 


z{t) < In - (7g - I^g) + 

I 7g V ^ 

y{ri) < In f(7 G-/xg) + 


(^NS 


T 


use^e ^dt = y.Ns'T 


x{C) < In fandx{^) > In f _ 

\ <^s J \ J 


Furthermore, 


(70) 


(71) 


(72) 


(73) 


J z{t)dt = j 

Using 

0 < (7g - Pg)'^ + ln(l - 7g) = 

< 


iG- l^G- 

J^G 


= -ln(l -7 g)- 


7G z(t) I aVd) 

0 L-^G 
'0 L^G 


dt 


dt = T 


J^G 


(74) 


then there exist Ci > 0, C 2 > 0 such that 

ln(l - t]g) 


Cl + C2 — (7g — Pg) + 

and 
Kg 


T 
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> (2- 


Thus, 


zir) > In <{ ^Ci } , y{v) > In 


{■ 


(ZnS 


■C2 


So, keeping in mind that by assumptions of Theorem [2l one has 


(7G-MG) + ihi^>0. 


x{t) < x{§) + \x{t)\dt 
< In 


Us 


< In 


y{t) < yi^) + lo \y{t)\dt 
I O'NS 

z{t) < z{t) + £ \z{t)\dt 
Kg 


+ {us + ys + 'ys)T + 

ln(l - yc) 


{iG — yc) + 


r 


:= 


• yu 


< In 


7g 


(7 g — yc) + 


ln(l - rjc) 


T 


+ 2-fGT + ln(l - yc) ■= Zu, 


x{t) > x{^) — \x{t)\dt 


> In 


yNs 

Us 


{us + ys + 75 )^ - - 'yNsTe^^^'^ 


y{t) > y{v)-£\y{t)\dt 

C2 


> In 


— 2y]s[sT 


_ O'NS _ 

z{t) > z{t) - \z{t)\dt 
'Kg 


> In 
therefore we obtain. 


7g 


Cl (■ - ‘^iG'r - ln(l - Vg) 


■ = xi, 


■= yu 


■= Zu 


max |a;(t)| < max||a;„|, IxiH := M^,, 

0<t<T n I IJ 

m^ \y{t)\ < max{|i/„|, \yi\} := My, 

0<t<T 

m^ k(^)| < niax{|^„|, \zi\} := 


Ma,, My and M^ are independent of /3. 

Now let us consider the algebraic equations 


1 

r 


7g 

Kt 


(e"" + e^) - /SyATge^e 


7s - czs - jUs + - ln(l - 

T 

\ 


—y-Ns + 


r 


(7g - fiG) + - ln(l - 7g) - - 


X Jo [Kg 


7G z : y 

-r^e +aNse^ 


Z{T)■^ 

dt = 0, 
use^e~'^dt = 0 , 
dt = 0 


(76) 


(76) 


(77) 


(78) 


for {x,y,z) G M^, where /3 G [0, 1] is a parameter. By carrying out similar arguments as in 
system dSZD, one can show that any solution {x*,y*,z*) of (175]) with (3 G [0, 1] satishes 


h<x* < Li, l 2 <y*< L 2 , h<z* <L 


3- 


(79) 
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Taking Mq = + My + + Mk where M^iO is taken sufficiently large such that 

Mk > |/i| + \Li\ + I/ 2 I + IT 2 I + I/ 3 I + l-hgl, we define Vt = [{x,y,zY' G X : \\{x,y,z)\\ < Mq}, 
then n verifies the requirement (1) of The Continuation Theorem (Gaines and Mahwin 
(1977) [27|, Page 40). When {x, y, z) G Oil fl KerL = dVt fl {x, y, z) is a constant vector 
in R^ with ||(x, y, ^)|| = Mq. Then from fl79l) and the definition of Mq, one has 



where, 


^( 1 ) = 1,^-1 £ + + + 
+ilii(l - i7sTO(e-'<^>)), 

T 

= —yNs + -( ujse^e~ydt, 

^ Jo r 

^( 3 ) = [ aNse^dt+-\n{l-r]G), 

r Jo Kg t Jo t 


(81) 

that is, the first part of (2) of The Continuation Theorem (Gaines and Mahwin (1977) [27], 
Page 40) is valid. 

To compute the Brouwer degree, let us consider the homotopy 


dt 


Hyiix,y,z)'^)= (3JQN{{x,y,z)'^) + {l-(3)V{{x,y,z)'^), /3 G [0, 1], 


where 


75 - wg - /i5 + - In (1 - yswiK^'"^)) - + e^) ^ 

1\T 


V{{x,y,zf) = 


V 


-Uns + ixse^e ^ 

(7g - /Wg) + - hi(l - 7 g) - 

T Kg 




y 




Yiom. fl78|l . it follows that 0 ^ Hp{dVLr\KerL) for j3 G [0, 1]. Moreover, since + ^^Ns) ^ 

KtKg 

0, the algebraic equation V{{x,y, z)'^) = 0 has a unique solution (e**, e^*, e^*)^ G R^. We 
compute the Brouwer degree {deg{-, •, •)) by using the invariance property of homotopy [?], 
one has 

deg{JQN,^l D KerL,0) = deg{V,^l D KerL,0) 

= E sign{Jv{p)) 


pey-i{o) 




( 

= sign 

det 




V 

= sign 

iGlS 

KtKg 

= -1^0. 



75 


Kt 

Pns 

0 


Kt 

—PNS 


-aNse^ 
(Pns + ^s)K 


0 

0 

7g^.* 


Kr 


\ 

/ 


, since ujse = Uns^ 


y 


(82) 
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By now, we have prove that f2 verihes all requirements of The Continuation Theorem (Gaines 
and Mahwin (1977) [2Z], Page 40), then Lx = Nx has at least one solution in DomL fl hi, 
i.e. system (15^ has at least one solution in HomLnH, say {x*(t), y*(t), z*. Set Tg(t) = 
exp(a;*(t)), T^si^) — 6xp(l/*(^)) G*it) = exp{z*{t)), then {Tg{t),T^g{t),G*{t))'^ is a positive 
and r-periodic solution of system ([I]) — ([2]). This completes the proof. 


Appendix D: Particular values of (^G r 


With respect to relations flT^ and flTTI) . we set 

1 


= 7^(7s-(/^ s + c^s + hvs)), 

^int 

_ (7s — (ls + 1^5 + Lns)) 

Kg f , ln(l -r]G)\ 

- 7g - hG H- 

7g V ^ 


’NS\ 


= 


1 f , ln(l - Vg) 

^ 7G-/iG +--- 

-L NS \ ^ 

ln(l - ijg) 




7g 1 


(83) 


One can note that aQ, and r* determined regions of stability/instability of forest and 
grassland solutions, with respect to ug, cns and r variations. 
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